APS/123-QED 



O 

(N 
IX, 



Origin of coexisting large Seebeck coefficient and metallic conductivity in the electron 

doped SrTiOg and KTaO s 

Hidetomo Usui, Shinsuke ShibataQ, and Kazuhiko Kuroki 

Department of Applied Physics and Chemistry, The University of Electro-Communication, Chofu, Tokyo 182-8585, Japan 

(Dated: February 10, 2010) 

We study the origin of the large Seebeck coefficient despite the metallic conductivity in the La- 
doped SrTiOa and Ba-doped KTaC>3. We calculate the band structure of SrTi03 and KTa03, 
from which the Seebeck coefficient is obtained using the Boltzmann's equation. We conclude that 
the multiplicity of the t2 g bands in these materials is one major origin of the good thermoelectric 
property in that when compared at a fixed total number of doped electrons, the Seebeck coefficient 
and thus the power factor are larger in multiple band systems than in single band ones because the 
number of doped electron bands per band is smaller in the former. We also find that the second 
nearest neighbor hopping integral, which generally has negative values in these materials and works 
destructively against the Seebeck effect, is nearly similar between KTaC>3 and SrTiOs despite the 
larger band width in the former. This can be another factor favorable for thermopower in the 
Ba-doped KTa0 3 . 
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I. INTRODUCTION 

The discovery of the large Seebeck coefficient in 
Na x Co02@ and the findings in cobaltates/cobaltites0- 
7] and rhodates [1, H| that followed have brought up an 
interesting possibility of finding good thermoelectric ma- 
terials that have relatively high (metallic) conductivity. 
These cobaltates and rhodates are materials where holes 
are doped into the d 6 configuration, namely the electron 
configuation where the t^g bands are fully filled. On 
the other hand, there is another class of tig transition 
metal oxides where relatively good thermoelectric prop- 
erties are obtained, namely the electron doped materials 
such as SrTiC>3 10]. When Sr is partially replaced by 
La in SrTi03, electrons are doped in the originally d 
configuration. This material exhibits large Seebeck co- 
efficient despite showing metallic conductivity, and the 
power factor, i.e., the Seebeck coefficient squared times 
the conductivity, is comparable to that of Bi2Te3. Quite 
recently, good thermoelectric pro perties have also been 
observed in Ba-doped KTaOajlJ]. This is another t2 g 
oxide, where electrons are doped into the originally d° 
configuration, but is a 5d system as compared to the 3<i 
in SrTiC>3. Here again, relatively large Seebeck coeffi- 
cient is observed despite the metallic conductivity. 

Theoretically, there have been several approaches that 
explain the large Seebeck coefficient in oxides. From 
the first principles band calculation studies, it has been 
pointed out that the narrowness of the band struc- 
ture in Na^Co02 and related rhodates is an important 
factor [T2I lill]. We have proposed that besides the width 
of the band, the shape of the band, which we call the 
"pudding-mold" type band, is important for the coexis- 
tence of the large Seebeck coefficient and the high con- 
ductivity in NazCoC^dH and related rhodates [H, [lj|. 
On the other hand, Koshibae et al. derived a formula 
for the Seebeck coefficient in the T(temperature) — > 
00 limit, and pointed out that the orbital degeneracy 



originates large entropy, leading to the large Seebeck 
coefficient 1 171 Il8|. 

In the present study, we propose that yet another 
mechanism, where the band multiplicity plays an impor- 
tant role, is at work in the electron doped t2 g materi- 
als. Namely, when there are multiple (nearly) equivalent 
bands at the Fermi level, and the number of doped elec- 
trons per band is fixed, the Seebeck coefficient is essen- 
tially the same regardless of the number of bands, while 
the conductivity increases with the multiplicity of the 
bands thus resulting in an enhanced power factor. In 
other words, when the total number of doped electrons 
itself is fixed, the Seebeck coefficient and thus the power 
factor is larger for multiple band systems because the 
Fermi energy stays low. We also examine the effect of 
the band shape, and show that the second nearest neigh- 
bor hopping integral, which generally has negative values 
in these materials and work destructively against the See- 
beck effect, is nearly similar between KTa03 and SrTi03 
despite the larger band width in the former. This can be 
another factor favorable for good thermoelectric proper- 
ties in the Ba-doped KTa03. 



II. FORMULATION 



A. Boltzmann's equation approach 



We first briefly summarize the Boltzmann's equation 
approach adopted in the present study [3, [l9j]. In this 
approach, the Seebeck coefficient is given as 



(1) 



where e(< 0) is the electron charge, T is the temperature, 
tensors K and are given by 



(a) SrTi0 3 



(b) KTa0 3 



K„ = r{k)v{k)v{k) 



9_m 

de 



(k) 



(e(k)-») n . (2) 



Here, e(k) is the band dispersion, v(k) — V^e(k) is the 

group velocity, r(fc) is the quasiparticle lifetime, /(e) is 
the Fermi distribution function, and /i is the chemical 
potential. Hereafter, we simply refer to (K n ) xx as K n , 
and S xx — (l/eT)(Ki/K ) (for diagonal K ) as S. Using 
Kf), conductivity can be given as a xx = e 2 K = a = 1/ p. 
As an input of the band structure in this calculation, 
we use the first principles calculation as described below. 
r(fc) will be taken as an (undetermined) constant in the 
present study, so that it cancels out in the Seebeck coef- 
ficient, while the conductivity and thus the power factor 
has to be normalized by a certain reference. 




FIG. 1: The band structure of (a) SrTi0 3 and (b) KTa0 3 In 
(a), the black dotted lines are the original LDA calculation 
while the solid red lines are the bands of the tight binding 
model obtained using the MLWFs. 



B. Band Calculation 

SrTiC>3 and KTa03 have cubic perovskite structures. 
We use the experimentally determined lattice constants 
in the band calculation, which are a = 3.90528A for 
SrTi0 3 [Io| and a = 3.9883A for KTa0 3 [E|]. For 
SrTiC>3, we have obtained the band structure using the 
Quantum-ESPRESSO package[22j. In order to obtain 
a tight binding model on which we can perform vari- 
ous analysis, we construct the maximally localized Wan- 
nier functions (MLWFs) [23| for the outer energy win- 
dow OeV < efe — Ep < 5eV and for the inner windows 
OeV < efc — Ep < 4eV, where ek is the eigenenergy of 
the Bloch states and Ep is the Fermi energy. These ML- 
WFs, centered at Ti sites in the unit cell, have three 
orbital symmetries (orbital l:d xy , 2:d yz , 3:d zx ). With 
these effective hoppings and on-site energies, the tight- 
binding Hamiltonian is obtained, and finally the See- 
beck coefficient is calculated using eq.([T|). For KTa03, 
we have obtained the band structure using the WIEN2K 
package [24) . The Seebeck coefficient is calculated using 
the BoltzTraP code[H]. 



III. CALCULATION RESULTS OF THE 
SEEBECK COEFFICIENT 

In this section, we present the band calculation results 
and the calculation results of the Seebeck coefficient. 
The calculated band structures of SrTiC>3 and KTa03 
are shown in FigfT] In both materials, there are three 
t2 g bands right above the Fermi level, and for SrTiC>3, 
the band structure of the three band tight binding model 
is superposed to the original first principles band. The 
band structure of the two materials look similar, but the 
band width is wider for KTa03 due to the widely spread 
nature of the 5d orbitals. 
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FIG. 2: The calculated Seebeck coefficient for (a) SrTiOs and 
(b) KTaC>3 plotted as functions of temperature for various 
doping rate x. 
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The calculated Seebeck coefficient for the two materi- 
als is shown in Fig(3]against the temperature at x = 0.05 
and x = 0.1 for SrTi0 3 , and x = 0.009 for KTa0 3 . 
We have chosen these x to make comparison with the 
experiments [ToL ITD] . Here we take a rigid band approach, 
and assume that the hole concentration nh is equal to the 
La (SrTi0 3 ) and Ba (KTa0 3 ) content. 

For SrTi0 3 . The Seebeck coefficient at 300K is S(x = 
0.05) = -87/xV/K and S(x = 0.1) = -58^V/K. Ex- 
perimentally, the Seebeck coefficient at 300K is S(x = 
0.05) = -147^V/K and S(x = 0.1) = -88.7/iV/K[10]. 
Thus the calculation result is somewhat reduced from 
the experimental result. The reason for this is proba- 
bly due to the reduction of the band width due to the 
strong correlation effect of the 3d orbitals. In fact, it has 
been known from the comparison between band calcula- 
tions and the angle resolved photoemission studies that 
the band width of the 3d electron materials is generally 
reduced by a factor of about two, and in fact taking this 
effect into account reproduces the experimental results of 
Na x Co02 well [lij]- If we calculate the Seebeck coefficient 
at 300K by reducing the band width by 50% from the 
bare LDA result, we get S(x = 0.05) = -149//V/K and 
S(x = 0.1) = -103^V/K(FigH;a)), which are in fact 
fairly close to the experimental values. 

As for KTa0 3 , the calculation of the Seebeck coeffi- 
cient at x = 0.009 gives S(300K) = -160/zV/K. This is 
roughly in agreement with the experimental result which 
is about — 200/zV/K[ll|. A reason why the bare LDA 
band structure gives good agreement with the experi- 
ments is because KTa0 3 is a 5d system, where the elec- 
tron correlation effects are expected to be small com- 
pared to 3d systems like SrTi0 3 . In fact, for a number 
of rhodates, i.e. Ad systems, the Seebeck coefficient cal- 
culated from the bare LDA band structure gives fairly 
good agreement with the experiments [lH. [TBI. Il6j . 



IV. EFFECT OF THE BAND MULTIPLICITY 

Having found that the experimentally observed See- 
beck coefficient is roughly reproduced within the first 
principles band calculation + the Boltzmann's equation 
approach (with some additional consideration of band 
narrowing), we now explain why the Seebeck coefficient 
is large in these materials despite the relatively large con- 
ductivity. In other words, we seek for the origin of the 
large power factor S 2 a. 

In the three orbital model, the Seebeck coefficient S xx 
is given as 



1 Kf xy + K dyz + K dzx 
eT K dxy + K$ yz + K dzx ' 



(3) 



where -K"„ stands for K n of the — x,y,z) or- 

bital. From eq.©, the group velocity u* J is the im- 
portant factor in K n . v dxy is equal to v because 



d{e xy )/dx is equal to d(e zx )/dx, so that K dxy 



K 



dzx 



Also, K dyz ~ because v dyz is very small. So the See- 
beck coefficient is 

1 OK dxy 1 K dxy 

eT 2K dxy ~ eT K dxy ~ xx ' 1 ' 

Namely, the total Seebeck coefficient is equal to the See- 
beck coefficient of the d xy single orbital system. On the 
other hand, the conductivity is 

a = e 2 {K dxy + K dyz + K dzx ) ~ 2e 2 K dxy = 2a dxy . (5) 

Therefore the power factor is 

I',, = oS 2 xx 2a dxy (S dx x y ) 2 = 2P dxy . (6) 

The left hand side here is the power factor of the three 
orbital system, while P dxy in the right hand side is that 
of the d xy single orbital system. Thus the multiplicity of 
the orbitals is advantageous for large power factor. Note 
that the comparison here between the three orbital and 
the single orbital systems is given for the same number 
of electrons per band. If we present this relation between 
the three and one orbital systems using the doping con- 
centration x, it should be giv6n as S X x (3x) — ^xx^ (*^0 
and P xx (3x) = 2P dxy (x). 

In FigUH we show the actual calculation result of the 
Seebeck coefficient and the power factor (normalized at 
x = 0.13 of the single orbital model) of the ti g three or- 
bital model of SrTi0 3 and a single orbital model where 
only the d xy orbital is considered. It can be seen that 
the above relation is indeed satisfied. It is also worth 
noting that the doping dependence of the power factor is 
in striking a gree ment with the experimental observation 
(Fig. 3 in ref.[10|). From this figure, we can see that for a 
fixed doping concentration, both the Seebeck coefficient 
and the power factor is larger for multiorbital systems 
than in single orbital ones. This can intuitively be un- 
derstood from FigJ3Jc), namely, the Fermi level tends to 
stay lower for systems with multiple bands for a fixed 
number of doped electrons, and lower Fermi level results 
in a large Seebeck coeffiecient, while the large number of 
electrons (due to the multiplicity of the bands) enhances 
the conductivity [lfj. The present result suggests that 
the band multiplicity is at least one of the main reasons 
why the Seebeck coefficient is large despite the metallic 
conductivity. The orbital degeneracy has been consid- 
ered as a factor to obtain good thermoelectric properties 
in the context of entropy [TJ [3, but we stress here that 
the present mechanism provides another way where the 
band multiplicity can play an important role[l6|. 



V. EFFECT OF THE BAND SHAPE 

In the present materials, the density of states per 
band is not so large around the Fermi level. This can 
roughly be understood in terms of the tight binding 
model. Namely, the tight binding model on a square 
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FIG. 3: (a) The Seebeck coefficient and (b) the normalized 
power factor of the single orbital model (dashed green) and 
the three orbital models (solid red) as functions of the doping 
rate x at 300K. (c) Schematic figure of how the Fermi energy 
differs between single and multiband sysmtems. 



lattice has electron hole symmetry when only the near- 
est neighbor hopping t\ is considered. The introduction 
of the second nearest neighbor hopping t 2 breaks this 
electron- hole symmetry, and for t 2g systems, this hopping 
integral usually has a negative sign when writing down 
the Hamiltonian in the form H = 23« *y c l c j- When 
t 2 is negative, the density of states tends to be large in 
the upper half of the band and small in the lower half. 
In this sense, the effect of the so-called "pudding mold 
type" band [14]], where a flat portion of the band has to be 
present near the Fermi level, is not relevant to the present 



TABLE I: ti, t2 and |t2/£i| obtained from constructing max- 
imally localized Wannier orbitals. In performing the band 
structure calculations, we have used the experimentally de- 
termined lattice parameters taken from the cited references. 
For KTa0 3 , we have obtained ti and ti by fitting the WIEN2k 
band structure with a tight binding model. 





ii(eV) 


ta(eV) 


l*a/*il 


PbTi0 3 [26] 


-0.23 


-0.073 


0.31 


BaTi0 3 [27j 


-0.25 


-0.066 


0.26 


SrTiO 3 [20j 


-0.28 


-0.078 


0.28 


BaZr0 3 [28j 


-0.40 


-0.081 


0.20 


NaNb0 3 [29] 


-0.45 


-0.091 


0.20 


KTa0 3 [21J 


-0.52 


-0.094 


0.18 


BaMnO 3 [30j 


-0.17 


-0.067 


0.41 



electron doped materials. This can in fact be seen as fol- 
lows. Since we have found that S xx ~ S x % y (assuming 
same electron number per band) in the preceding sec- 
tion, we concentrate here on the d xy single orbital model 
of SrTiC>3. In this model, the nearest and second nearest 
neighbor hoppings of the MLWF tight-binding Hamilto- 
nian are t\ = — 0.28eV and t 2 = — 0.078eV. To see how 
t 2 affects the Seebeck coefficient, we vary t 2 while fix- 
ing t\ = — 0.28eV, and calculate the Seebeck coefficient 
at 300K as shown FigQJa). It is found that the smaller 
1*2 1 is, the larger the Seebeck coefficient. This is because 
the lower part of the band (where the Fermi level ex- 
ists) become less dispersive as \t 2 \ is decreased when t 2 
is negative. This can be seen in the calculation of the 
density of states (DOS) given in Figj4jb), namely, the 
DOS at the band bottom for t 2 — OeV is about twice 
larger than for t 2 = — 0.13eV. Thus the negative value 
of t 2 (i.e., the band shape) in SrTi03 is not faverable 
for thermopower, and the good thermoelectric properties 
seem to come mainly from the multiplicity of the bands. 

We have also evaluated t\ and t 2 for KTa03 from the 
obtained band structure as listed in table U together with 
some related materials. Although t\ is much larger com- 
pared to that in SrTiOa as expected from the 5d nature, 
t 2 is not much changed, and the ratio \t 2 /t\ | is the small- 
est among the materials considered here. In fact, \t 2 /ti\ 
is also small in Zr and Nb compounds, namely 4d systems 
with small number of electrons. So it seems that the ra- 
tio |ta/ii| tends to be small for large principle quantum 
number. This trend can be considered as another factor 
working favorable for the thermopower in KTa03 despite 
the wide band width. 



VI. CONCLUSION 

To conclude, we have studied the origin of the large 
Seebeck coefficient in SrTi0 3 and KTa0 3 . In SrTi0 3 , 
from the first principles band calculation results, a tight- 
binding model is obtained via the maximally localized 
Wannier orbitals, and the Seebeck coefficient is calcu- 
lated using the tight-binding model. In KTa03, from 
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FIG. 4: (a) The Seebeck coefficient of the single orbital model 
for x = 0.017, ti = -0.28eV and T = 300K plotted as 
functions of t^. (b) The density of states at t\ = — 0.28eV, 
t2 = OeV (solid red) and ti = — 0.13eV (dashed green). 



the first principles band calculation results, the Seebeck 
coefficient is calculated using the BoltzTraP code. In 
both materials, the large Seebeck coefficient despite the 
relatively large conductivity is largely due to the mul- 
tiplicity of the bands, i.e., essentially the same value of 
the Seebeck coefficient is obtained for the same number 
of electrons per band, so that when the total number of 
doped electrons itself is the same, the Seebeck coefficient 
and thus the power factor are larger for multiple band 
systems. Also, we have examined the effect of the band 
shape. Although the negative t% value is not favorable 
for the electron doped thermoelectric materials, 4d and 
5d systems such as KTa03 tend to have similar ti val- 
ues as in 3d systems despite the wide band width, and 
this can be another factor that is advantageous for good 
thermoelectric properties. 
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